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The exceptional capability of plasmonic structures to confine light into deep subwavelength volumes has 
fashioned rapid expansion of interest from both fundamental and applicative perspectives. Surface plasmon 
nanophotonics enables to investigate light - matter interaction in deep nanoscale and harness 
electromagnetic and quantum properties of materials, thus opening pathways for tremendous potential 
applications. However, imaging optical plasmonic waves on a single nanometer scale is yet a substantial 
challenge mainly due to size and energy considerations. Here, for the first time, we use Kelvin Probe Force 
Microscopy (KPFM) under optical illumination to image and characterize plasmonic modes. We 
experimentally demonstrate unprecedented spatial resolution and measurement sensitivity both on the 
order of a single nanometer. By comparing experimentally obtained images with theoretical calculation 
results, we show that KPFM maps may provide valuable information on the phase of the optical near field. 
Additionally, we propose a theoretical model for the relation between surface plasmons and the material 
workfunction measured by KPFM. Our findings provide the path for using KPFM for high resolution 
measurements of optical plasmons, prompting the scientific frontier towards quantum plasmonic imaging 
on submolecular scales. 

Surface Plasmon Polaritons (SPPs) are optically excited collective charge oscillations, which can propagate 
along metal - dielectric interface 1 . These electromagnetic waves are characterized by subdiffraction volumes 
and high intensities, properties which unlock the potential for huge variety of applications 2 " 6 . Nevertheless, 
since being bound to the surface and confined to deep subwavelength scales, high resolution imaging of nano- 
plasmonic modes using conventional optical microscopy is nearly impossible 7 " 9 . Subdiffraction optical imaging of 
surface plasmons can be achieved using near field techniques, as scanning near field optical microscope (SNOM) 
and its modifications that use metal tips 10 " 15 . In SNOM, a subwavelength optical aperture samples the average 
optical signal intensity as a function of position. However, spatial resolution and collection efficiency are both 
limited by the aperture size and by the distance between a sample and the tip. This is mainly because SNOM 
samples only the decaying or evanescent part of the optical signal outside the nanoscale object 16 . To date smallest 
commercially available SNOM tips achieve aperture diameters of 50 nm - 100 nm (Nanonics - Optometronic 
4000™ at http://www.nanonics.co.il/products/nsom-spm-systems/optometronic-4000.html). Phase resolved 
imaging of optical near field with resolution on the order of 10 nm - 20 nm was recently demonstrated using 
Heterodyne and pseudoheterodyne scattering type SNOM 17 (s-SNOM). Plasmonic imaging with nanometric 
spatial resolution can also be achieved using electron energy loss spectroscopy (EELS) 18-20 , with time resolved 
measurements are also demonstrated using laser excited photo cathodes in electric microscopy 21 ' 22 . Kelvin probe 
force microscopy (KPFM), also known as surface potential microscopy, is a variant of atomic force microscopy 
(AFM) which was first introduced in 1991 23 . KPFM enables to measure the work function difference between a 
scanning tip and a surface with submolecular resolution and submilivolt sensitivity 24 " 28 . When a scanning probe 
tip is electrically connected to a conductive sample, a contact potential difference (CPD) will arise due to the 
different workfunctions of the tip and the sample. KPFM measures local variations in the CPD by applying a 
voltage between the sample and the tapping AFM tip so that the electric field caused by the CPD and the resulting 
force on the tip are compensated. For a certain tip position, the compensating voltage represents the local contact 
potential difference (LCPD) and can be determined either with a feedback loop or by measuring the voltage 
dependence frequency shift 23,26 ' 29 " 31 . Recently, we used KPFM and space- time Fourier analysis to characterize a 
nanoplasmonic device with resolution of 10 nm 32 . However, no imaging mechanism of SPPs with KPFM has so 
far been presented, and achieving single nanometer resolution imaging of optical fields is an extraordinary 
practice. Here, we show that KPFM can image and characterize surface plasmon polaritons with resolution at 
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the single nanometer scale. A theoretical model is proposed to define 
the relation between CPD and SPPs and provide an in-depth 
understanding of the imaging mechanism. We characterize plasmo- 
nic nanoantenna and waveguiding structures achieving measure- 
ment sensitivity of only 2 nm. Additionally, we observe that the 
sign of the CPD can be used to obtain information on the electric 
field's directionality, towards full vector characterization of SPPs. 
Figure 1 illustrates the process of SPP characterization with scanning 
KPFM. The origin of nanoplasmonic sensing with KPFM lies in the 
fact that KPFM measures the contact potential difference (CPD) 
between the scanning probe tip and the surface 26 ' 33 . The CPD is 
closely related with the charge distribution on the surface, as well 
as with the surface plasmon polaritons 34 " 40 . When nanometallic 
structures are illuminated with optical radiation, coherent oscilla- 
tions of free electrons are excited on the surface. These highly con- 
fined optically induced charge density waves preserve the frequency 
of incident photons and strongly affect the work function in the 
vicinity of the device. Optically excited variations in the workfuction 
of CdSe nanowires 41 and 2D array of subwavelength apertures 42 was 
recently characterized by EFM and KPFM. 

Results 

According to the Lorentz model, the motion of electrons in metals is 
governed by a second order partial differential equation 



d x ^dx ~ 



-eE 



(1) 



Here m and e is the electron mass and charge, respectively; x is the 
electron displacement from the nucleus T = 1/ T is the damping rate, 
Q is the spring constant and E is the electric field. Assuming time 
harmonic behavior for the electric field and displacement, neglecting 
Q in the free electron region yields the Drude polarizability: 

Ne 2 1 , . 

-E(x,t) (2) 
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where co = 2nf is the optical angular frequency of the electric field. 
The in plane relative permittivity (e| = s r + ie,-) can be extracted from 
Equation 2, using s = 1+ P/E. The real and imaginary part of 8\\ are 
attributed intraband and interband transitions, respectively. 
Extended Drude model which considers the restoration force (Q) 
is a preferred analytical description at the bound electron region. 
This model achieves excellent agreement for Au within 500 nm- 
1000 nm wavelength range 43 . Further improvement to the extended 
Drude model can be obtained by adding terms which carefully 



describe the critical points (CP) 44 ' 45 . Near the plasma frequency, 
co~co p = \/Ne 2 /ms 0 , the relative permittivity is very small 
(s(cop) ->0). Hence, E(x, t) becomes a pure depolarization field with 
quanta of these charge oscillations are referred to as plasmons. Using 
the relation between polarization and charge density p pol = — VP, 
the surface electric potential is calculated by integrating the charge 
density function over the computational surface 
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The charge density function may also be obtained by solving the 
Fock equations using an exchange potential operator 46,47 , a more 
accurate solution but computationally expensive when obtained 
numerically. The experimental setup used in this work is depicted 
in Figure 2a. The sample is illuminated from above by unfocused He- 
Ne laser at wavelength of 633 nm (<x> = 2.98T0 15 r/s) at incident 
angle of 45°. The light beam is linearly polarized in parallel with 
the dipole orientation, with the sample is appropriately oriented to 
avoid the cantilever from blocking the light path to the devices. Since 
co>r~10 13 , the in plane relative permittivity can be written as 

-.2 



co 2 / co 2 + icOp I co z t, with s B is the contribution of bound 



electrons and co p is the plasma frequency. In our measurements 
topography and KFPM images are recorded sequentially using the 
lift mode technique to minimize cross talk between the measure- 
ments 48 . To this end, we first acquire the surface topography of a 
single line in tapping mode and then retrace this topography over the 
same line at a set lift height from the surface to measure the KPFM 
signal. To map the workfunction of the sample, we apply both AC 
voltage (V A c) an d a DC voltage (V DC ) to the AFM tip. V A c generates 
oscillating electrical forces between the AFM tip and sample surface, 
and V DC nullifies the oscillating electrical forces that originated from 
CPD between tip and sample surface. In our study, the tip - sample 
distance (H) significantly influence the resolution and sensitivity of 
the measurements. For small values of H, the KPFM image is domi- 
nated by short range forces (~H~ 6 ) like Van der Waals and capillary 
interactions. However, for larger values of H, long range electrostatic 
and magnetic interactions (~H~ 2 ) govern the results. Therefore, it is 
critical to optimize the lift height for high quality imaging. To dem- 
onstrate unambiguously that we can measure and image plasmonic 
waves with nanoscopic resolution and sensitivity we chose to invest- 
igate hybrid devices of dipole nanoantennas integrated with metal - 
insulator - metal (MIM) plasmonic waveguides as shown in 
Figure 2b. These devices are of prime importance for large variety 



Nanoantenna Coupled to 
MIM Plasmonic Waveguide 




Figure 1 | Illustration of nanoplasmonic characterization with KPFM. Dipole nanoantenna is connected to MIM waveguide with nanometric gap 
width. The device is illuminated by a laser and the excited plasmons change the CPD between the tip and the surface, which measured by KPFM with 
nanometric resolution. 
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Figure 2 | Schematics description of our KPFM experimental process and fabricated devices, (a), Schematic of the measurement principle. First the 
topography is acquired in tapping mode, then the tip elevates to a set lift height, H to create the CPD map. (b) , Description of the analyzed nanoplasmonic 
devices. Dipole nanoantenna with arm length 'L A ' and gap 'g', are connected to MIM waveguide of length L w . The thickness and width of all metal 
strips are 12 nm. 



of applications 49,50 , and have several advantages for this study. First, 
the combination of nanoantenna and plasmonic waveguide enables 
to image both localized and propagating SPPs. Localized SPPs are 
optically excited at the nanoantenna and are coupled to propagating 
modes along the waveguide. Furthermore, by comparing CPD 
images with calculated vector fields we show that KPFM measure- 
ments may provide information on the vector optical near field. The 
devices where fabricated by electron beam lithography (EBL), ion 
beam sputtering (Ag, 10 nm) and lift off. By fabricating the struc- 
tures on the insulating layer of a Si substrate (2 |im thick Si0 2 ), we 
avoid coupling effects between the conducting Si and the metallic 
device. After the lift off the PMMA is completely removed, allowing 
contact between the AFM tip and the metallic devices. Figure 3a 
presents a scanning electron microscopy (SEM) image of the 




analyzed device, with geometrical dimensions of L w = 750 nm, L A 
= 100 nm, and g = 12 nm. The image was recorded at beam current 
of 0.4 nA and low accelerating voltage of 5 kV. Figure 3b shows an 
AFM topography image of the device in 3D, obtained using Si tip 
with 2 nm typical diameter. Subsequently, KPFM characterization of 
the same device was performed in lift mode while the device is 
illuminated by a He-Ne Laser. The KPFM measurements were per- 
formed for tip - sample distances ranging from 10 nm-250 nm, with 
the stage remains spatially fixed and the tip is scanned. We observe 
that for small values of H, short range interactions strongly affect the 
KPFM image. For larger tip-sample distances, the results are gov- 
erned by the long range interactions. The range of distances for 
which we obtained the best resolution was between 30 nm-50 nm. 
Figure 3c shows the KPFM image recorded for tip-sample distance of 
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Figure 3 | SEM, AFM KPFM and numerical analysis of MIM nanoplasmonic device, (a), High resolution SEM image of the fabricated and analyzed 
device (b), 3D AFM topography image of the device. The image was recorded with Si super sharp tip with 2 nm aperture, (c), CPD image of the device 
under optical illumination with He-Ne laser at frequency of 474 THz, recorded with the same tip as (b). KPFM Signal scale bar: — 1.7 V-+ 1.7 V. 
(d), 3D Numerical calculation (FEM) results of the of the analyzed device under optical illumination of He-Ne laser, similar to the experimental setup 
used to obtain the results in (c). Scale bar: 50 nm. 
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Figure 4 | Energy transport and propagation properties of dipole nanoantenna coupled to MIM nanoplamonic waveguide, (a), 3D Schematics of the 
analyzed devices (b), Numerically calculated near-field image showing Re(E z ) = I E z l cos(4> z ). (c), Numerically calculated near-field image showing Re(E x ) 
= lE x lcos(4> x ). (d), Numerically calculated near-field image showing Re(E y ) = lE y lcos(4> y ) } with the inset shows E y for normalized field scale of ± 1/4 
(instead of ±1). 



H = 35 nm, which exhibits clear nanometric resolution. We observe 
two dipolar-like modes on the antenna segments, with strong fields at 
the antenna ends and at the gap, as expected for a dipole plasmonic 
nanoantenna in the visible 32 . The localized fields at the nanoantenna 
gap extend along the MIM waveguide, periodically changing their 
polarity. The modal behavior of the propagating SPPs along the 
waveguide is observed in the KPFM image (Fig. 3c). Four periods 
of modes appear along the waveguide, represented by peaks of the 
KPFM signal, as the measured 32 SPP wavelength is 240 nm. The 
shape of the waveguide, as well as the field decay into its metallic 
arms are very well captured. This provides direct experimental 
imaging of a propagating surface plasmon polaritons, obtained for 
the first time for channel waveguides with 12 nmX 10 nm gap cross- 
ection. Strong localized fields are observed also on the exterior metal 
- insulator interfaces of the waveguides, as predicted by the theory of 
SPPs 1 . We observe that the measured KPFM fields on the ends of the 
nanoantenna are inversely polarized, with the field on the left edge is 
positive +1.5 V and the field on the right edge is negative —1.5 V, 
corresponds to local optical near field direction of +z and — z, 
respectively (see Fig. 3d). This behavior is reproduced in numerical 
calculation results, presented at the optical frequency of 474 THz 
(633 nm). The numerical results are obtained using Ansys High 
Frequency Structure Simulator (HFSS V15 32 ' 51 ' 52 ) based on the finite 
element method (FEM). Generally, the electric field is described by a 
three dimensional vector E = (E x , E y , E z ), where each field compon- 
ent Ej is characterized by both magnitude |E f | and a phase (|)i 53 . 
Figure 3d presents the 3D electric field vector, |E|, with each arrow 
represents the direction of the local electric field vector. In spherical 
coordinates system, the angles 6, (p of each arrow shown in Figure 3d 
are the angles define the direction of the local near field vector. 
Figure 4 shows the numerically calculated scalar components of 
the electric fields on the device. Figure 4b shows the real part of 
the vertical near-field component Re(E z ) = |E z |cos((() z ) where |E Z | 
is the near-field amplitude and <\> z the phase. We observe high intens- 
ity at the nanoantenna ends and gap as well as periodical change in 
the field polarity, as expected for optical gap nanoantennas 17 . For the 



horizontal near field component E x , as shown in Figure 4c. For 
Re(E x ), we can see that the amplitude image exhibits a completely 
different pattern, featuring a highly localized, intense fields at the 
nanoantenna gap and ends with a spatially constant phase 17 . 
Similarly, for the E x near field component, there is no phase differ- 
ence between the two waveguide arms. Since E x has significantly 
higher amplitude than E z , the phase behavior inside the waveguide 
is dominated by E x . This result explains why no phase difference 
between the fields inside the waveguide is observed in the KPFM 
image (Fig. 3c). As expected, the E y component of the near field 
exhibits much lower enhancement, see Figure 4d. Interestingly, the 
phase structure of E y may explain the slight shift of the phase center 
of the KPFM signal (Fig. 3c) with respect to the topography. Figure 5a 
presents a high resolution AFM image of the fabricated device with 
gap of 4 nm (g = 4 nm, L A = 20 nm), which sets the lower limit for 
the gap width achievable via our EBL fabrication process. The KPFM 
signal of the corresponding topography is shown in Figure 5b. The 
device exhibits lower field enhancement, which may be attributed to 
the off resonance dimensions. Additionally, the gap width may 
approach ~1 nm at some points due to undesired roughness of 
the lithography process. This may lead to decreased field enhance- 
ment originated from quantum charge transfer 54 . Figure 5c shows the 
topography (blue) and KPFM (black) of the device plotted along the 
dipole axis. We observe that the KPFM signal of the different 
nanoantenna arms can be clearly resolved, while separated by only 
~4 nm gap. This provides first direct experimental evidence of 
KPFM capability to image surface plasmon modes with unpreced- 
ented resolution on the order of a single nanometer. This imaging 
resolution of KPFM holds potential for further improvements which 
can be obtained using sharper tips and process optimization 24 Now, 
we demonstrate another substantial vantage of KPFM for plasmonic 
characterization, by resolving the propagation properties of MIM 
waveguides with gap (g) difference of 2 nm. Figure 6a presents the 
configuration of the analyzed devices. Figure 6b shows the calculated 
electric field at frequency of 474 THz (k = 633 nm,) for devices with 
L w = 900 nm, g = 10 nm, 20 nm, 30 nm and 60 nm. We observe 
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inverse relation between the wavelength of the propagating SPPs, 
^ SPP , and the gap width. As the gap width decreases, the plasmons 
are confined in both dimensions and vice versa; a confinement which 
increases the propagation loss. This nonlinear dispersion behavior is 
evidenced by the difference number of SPP periods propagate within 
the waveguide. For g = 10 nm, seven periods of SPP modes are 
contained inside the waveguide, as for g = 60 nm only four periods 
are observed (see Fig 6b). The calculated electric field along the center 
of the waveguide (dashed line in Fig. 6a) is shown in Figure 6c, which 
further demonstrate the change in ^ SPP as the gap width increases. 
Figure 6d shows the KPFM signal map for g = 30 nm and tip - 
sample distance of H = 35 nm, selected for maximum resolution. 
This image has similar characteristics to Figure 3c with two main 
differences. First, only three periods of SPPs are contained within the 
waveguide (compared with four periods in Fig. 3c), this is due to the 
increase in gap width from 12 nm to 30 nm; a result which is in very 
good agreement with our calculations (Fig 6b). Second, lower field 
intensities are observed at the gap and on both nanoantenna ends, 
mainly caused by impedance mismatch between the dipole and the 
waveguide and by geometrical shift from resonance. The exact spatial 
frequency of the propagating SPPs is detected by analyzing the 
frequency content of the KPFM measurements along the waveguide 
via spatial Fourier transform (FFT). The black graph in Figure 6e 
presents the results of spatial FFT performed on the KPFM signal 
inside the waveguide (marked by line 1 in Fig. 6d), with the blue 
curve corresponds to spatial FFT performed outside the device (line 
2). Both measurements shown in Figure 6e have dominant field 
component at low frequencies, a result which is in good agreement 
with the general electrostatic character of KPFM and can be referred 
to the standing waves in the structure. However, we observe strong 
spatial frequency component in the field along the device, X SPP = 
248 nm. The amplitude of this component is 15dB higher than the 



corresponding component in the measurement described by the blue 
curve, providing solid evidence for the existence of propagating SPPs 
inside the waveguide. We used this method to characterize the pro- 
pagation properties is a series of fabricated nanoplasmonic devices 
with gap width varying between 10 nm-30 nm in steps of only 2 nm 
(i.e. g = 10 nm, 12 nm,. . .30 nm), L A = 100 nm and L w = 3000 nm. 
For each KPFM measurement, we found X SPP by performing spatial 
FFT and measuring the plasmon - wavelength at the peak of the curve 
(see Fig. 6e).We use long waveguides (L w > ^spp) in order to min- 
imize the standing waves in the device and improve the sensitivity. 
Figure 7a shows the experimentally obtained dispersion curve (blue 
graph), with the calculation results presented in the black curve. Real 
time image of the characterization process via KPFM under laser 
illumination if presented in Figure 7b. We observe different experi- 
mentally obtained plasmon propagation wavelength (^spp) for wave- 
guides with gap difference of 2 nm, which confirm that KPFM can 
resolve SPPs with 2 nm resolution, currently limited by the aperture 
dimensions of the scanning probe. The maximum deviation between 
the calculated and experimental dispersion curves is 5 = 2.5 nm, less 
than 1% from the mean ^ S pp- 

Discussion 

In Summary, we introduced KPFM as an exceptional method for 
imaging and characterization of plasmonic modes at optical frequen- 
cies. By investigating channel waveguides coupled to nanoantennas, 
we characterized both surface plasmons and surface plasmon polar - 
itons at the nanoscale. We characterized a variety of nanoplasmonic 
devices, with experimentally demonstrated spatial resolution and 
sensitivity both on the order of a single nanometer. The KPFM 
images include information on the amplitude and phase of the opti- 
cal near field, and are of good agreement with full wave numerical 
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Figure 6 | Spatial frequency analysis of the nanoplasmonic devices using KPFM, spatial FFT and numerical calculations, (a), Device configuration, 
(b), Numerically calculated magnitude of the optical near field, IEI devices with g = 10 nm, 20 nm, 30 nm and 60 nm, under optical illumination. 
Decrease in the spatial frequency (higher ^ SPP ) is observed when g is increased, (c), Calculated I El at 474 THz, plotted along the center of the MIM 
waveguides (dashed line in a) (d), CPD image of the device under optical illumination recorded at H = 35 nm. KPFM Signal scale bar: ± 1.2 V (e), Spatial 
FFT results oftheKPFM signal inside (black graph, calculated along Line 1 in (d)) and outside the waveguide (Line2 in (d)). Scale bars: 100 nmin (a) and 
50 nm in (d). 
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Figure 7 | SPP Characterization with 2 nm resolution, (a), ^ SPP as a function of the gap width, g, obtained experimentally (blue dots) and theoretically 
(continuous black curve), (b), Real time image of the characterization process, captured from the AFM system. Each square includes nanoantennas 
and MIM waveguide with different gap width. The laser polarization is aligned with the dipole orientation. 
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calculations. Our findings open the way for using KPFM under 
optical illumination for direct imaging of SPPs with nanometric 
spatial resolution and sensitivity, establishing KPFM as a preferred 
candidate for experimental characterization of SPP nano devices. It 
will now be possible for broad community of KPFM users to invest- 
igate optically induced charge density waves with dimensions com- 
patible to state of the art future CMOS technology nodes. 

Methods 

AFM and KPFM measurements. All measurements were performed at room 
temperature and free ambient conditions (no vacuum), using Dimension Icon AFM 
system with NanoScope V controller (Bruker®). For both AFM and KPFM 
measurements, we used Nano World probes SSS-NCH, SuperSharpSilicon - Non- 
Contact/Tapping™ mode - High resonance frequency; with typical diameter of 2 nm, 
resonance frequency of 320 kHz and spring constant of 42 N/m. typically, voltages of 
2 V, ac capacitance frequencies of 880 MHz, lift heights of 30 nm-50 nm and line 
rates of 0. 1 KHz were employed. 

Numerical simulations. The numerical results are obtained by using the software 
package ANSYS HFSS™ VI 5, the industry- standard simulation tool for 3D full-wave 
electromagnetic field simulation. HFSS solve Maxwell's equations via the finite 
element method (FEM) using adaptive mesh refinement process for tailored accuracy 
requirements. The field's solutions are calculated with the metallic (Ag) plasmonic 
structures being deposited on a homogenous Si02 substrate. The nanoantenna is 
illuminated by optical sources at 474 THz (wavelength of 633 nm), which are 
modeled as focused Gaussian beams with 1 um characteristic diameter. The electric 
field is polarized in parallel with the dipole direction, as the wave vector K is 
perpendicular. A selectively dense meshing is assigned in the metallic and 
waveguiding regions, with maximum cell size of 1 nm and 750,000 FEM tetrahedral 
cells. To provide maximum accuracy, the model is terminated as following: the 
interface with free space is bounded by perfectly matched layer (PML) absorbing 
boundary conditions (ABC), while the metallic and Si02 termination are done via 
layered impedance (LI) ABC. The minimum number of adaptive meshing iterations 
was set to 12, with convergence condition of 1% maximum energy variance between 
adjacent iterations. 

Devices Fabrication. Si02/Si sample was spin-coated with poly (methyl 
methacrylate) (PMMA 950 A2) electron-beam resist providing thickness of 100 nm. 
The samples coated with PMMA were subsequently baked for 120 s on a hotplate at 
180C. The desired pattern was exposed in the PMMA layer using a CRESTEC 
CABLE-9000C high-resolution electron-beam lithography system using different 
doses to control line and gap width. Then the samples were developed for 90 sec using 
methyl isobutyl ketone (MIBK), and rinsed with IPA. The samples were subsequently 
exposed to Ar plasma to etch 10 nm in order to remove leftovers from the pattern, 
sputtered using BESTEC 2" DC magnetron to deposit 3 nm Cr, and 13 nm Au, then 
immersed in 180 Khz ultrasonic bath with NMP for 3 h for resist liftoff. 
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